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Abstract 

Purpose Until recently, life cycle assessments (LCAs) have 
only addressed the direct greenhouse gas emissions along a 
process chain, but ignored the C0 2 emissions of land- 
use. However, for agricultural products, these emissions 
can be substantial. Here, we present a new methodology 
for including the implications of land occupation for C0 2 
emissions to realistically reflect the consequences of 
consumers’ decisions. 

Method In principle, one can distinguish five different 
approaches of addressing the C0 2 consequences of land 
occupation: (1) assuming constant land cover, (2) land-use 
change related to additional production of the product under 
consideration, (3) historic land-use change, assuming histor¬ 
ical relations between existing area and area expansion (4) 
land-use change related to less production of the product 
under consideration (“missed potential carbon sink” of land 
occupation), and (5) an approach of integrating land con¬ 
version emissions and delayed uptake due to land occupa¬ 
tion. Approach (4) is presented in this paper, using LCA 
data on land occupation, and carbon dynamics from the 
IMAGE model. Typically, if less production occurs, agri¬ 
cultural land will be abandoned, leading to a carbon sink 
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when vegetation is regrowing. This carbon sink, which does 
not occur if the product would still be consumed, is thus 
attributed to the product as “missed potential carbon sink”, 
to reflect the C0 2 implications of land occupations. 

Results We analyze the missed potential carbon sink by 
relating land occupation data from LCA studies to the po¬ 
tential carbon sink as calculated by an Integrated Global 
Assessment Model and its process-based, spatially explicit 
carbon cycle model. Thereby, we account for regional dif¬ 
ferences, heterogeneity in land-use, and different time hori¬ 
zons. Example calculations for several livestock products 
show that the C0 2 consequences of land occupation can be 
in the same order of magnitude as the other process related 
greenhouse gas emissions of the LCA, and depend largely 
on the production system. The highest C0 2 implications of 
land occupation are calculated for beef and lamb, with beef 
production in Brazil having a missed potential carbon sink 
more than twice as high as the other GHG emissions. 
Conclusions Given the significant contribution of land oc¬ 
cupation to the total GHG balance of agricultural products, 
they need to be included in life cycle assessments in a 
realistic way. The new methodology presented here reflects 
the consequences of producing or not producing a certain 
commodity, and thereby it is suited to inform consumers 
fully about the consequences of their choices. 

Keywords Agriculture • Carbon sink • Land occupation • 
land-use • LCA • Livestock 


1 Introduction 

Recently, the contribution of the livestock sector to green¬ 
house gas emissions and global warming has attracted a lot of 
attention (Steinfeld et al. 2006), and its potential contribution 
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to climate policy has been analyzed (Stehfest et al. 2009). 
According to the FAO, emissions from the livestock sector 
amount to 18 % of global greenhouse gas emissions in the past 
decade. The total GHG effect of the livestock sector is com¬ 
posed of the emissions during the production process, and the 
emissions related to land-use change due to the expansion of 
agricultural area and—as we will discuss here—also due to 
the use of agricultural land. While emissions along the pro¬ 
duction chain can be estimated quite accurately, the emissions 
from land-use are highly uncertain. 

Most of these estimates are based on so-called life cycle 
assessment (LCA), which is the state-of-the art methodolo¬ 
gy to calculate a product’s impact on the environment. The 
LCA methodology has been standardized by the Internation¬ 
al Organization for Standardization (ISO 14040 and 14044, 
see (ISO 2006)), and is further specified for the livestock 
sector (Curran 1993; Hendrickson et al. 1998; Guinee et al. 
2002; Gerber et al. 2010). All inputs and outputs along the 
production chain are compiled in a so-called life cycle inven¬ 
tory, and then combined to impact categories like total GHG 
emissions or global warming. The hitherto existing LCA 
results for livestock products for beef vary from 11 to 36 kg 
C0 2 -eq/kg (Ogino et al. 2007; Casey and Holden 2006; 
Williams et al. 2006; Blonk et al. 2008; Hirschfeld et al. 
2008) and for milk mostly from 0.9 to 1.5 kg C0 2 -eq/kg 
(Haas et al. 2001; Cederberg and Flysjo 2004; Casey and 
Holden 2005; Forster et al. 2006; Thomassen et al. 2007; 
Hirschfeld et al. 2008), and up to 7.5 kg C0 2 -eq/kg for sub- 
Saharan Africa (Gerber et al. 2010). 

Until recently, LCAs have not addressed the various 
effects of land-use on the climate system via modified fluxes 
of C0 2 and non-C0 2 gases, modified albedo, and modified 
evapotranspiration. Among these, modified fluxes of C0 2 
are the most relevant, and will be the focus of this study. 
While conventional LCAs do include direct emissions along 
a specific production chain in steady-state-conditions, it was 
not common practice to include carbon emissions or carbon 
sinks related to land-use change and land occupation (Mila i 
Canals et al. 2007). More recent papers either recognized 
land-use change as a relevant issue, but did not integrate it in 
the LCA results (Hirschfeld et al. 2008) or integrated his¬ 
toric changes in land-use into LCA (e.g., Gerber et al. 
(2010) for FAO). However, with global emissions from 
land-use change amounting to about 20 % of total green¬ 
house gas emissions (Rogner et al. 2007), they obviously 
must be included in the LCA for all products which require 
land (e.g., crops, livestock, bio-energy). The methodological 
challenge of doing so is still not solved. For bio-energy 
crops, with their fast expanding production and purpose to 
reduce GHG emissions, the relevance of land occupation 
and the related (direct and indirect) land-use change emis¬ 
sions has been studied extensively in recent years (Search- 
inger et al. 2008; Fargione et al. 2008; Melillo et al. 2009; 


Plevin et al. 2010; Overmars et al. 2011). For food crops and 
livestock products, however, the development was much 
slower. However, also for these other commodities, land 
occupation and land-use change have a huge impact on global 
warming and should be included into LCAs, too. As we show 
in this paper, land occupation itself (regardless of real 
changes in land-use) affects global warming as it prevents 
natural vegetation from regrowth and thus from carbon up¬ 
take (the “missed potential carbon sink”). In very general 
terms, it had been suggested earlier that the occupation of 
land might be understood as the prevention of regeneration 
(Doka et al. 2002). A recent paper already suggests a general 
methodology for including carbon implications from land 
occupation and land conversion into LCAs (Muller-Wenk 
and Brandao 2010). The differences to the method presented 
here will be discussed within this paper. 

In principle, the following five approaches and their 
inherent assumptions can be distinguished (Table 1): (1) 
current average of production, assuming a static system 
(no land-use change emissions, most conventional LCA 
studies until now). (2) Additional production of the product 
under consideration requiring additional land, which needs 
to be converted (mostly applied for biofuels, equally spread¬ 
ing conversion emissions over 30 years (Searchinger et al. 
2008), (3) average conversion emissions across a certain his¬ 
torical period, with mostly increased production, thereby in¬ 
cluding both the use of existing land and expansion (this 
method has e.g. been used by the FAO, see (Steinfeld et al. 
2006), and integrated in LCA by the FAO (Gerber et al. 2010), 
results strongly depend on increase of production in the respec¬ 
tive period). A recent paper (Ponsioen and Blonk 2011) sug¬ 
gests such a method specifically per country and crop type, 
allocating recent loss of natural land between forestry and 
agriculture, and between different agricultural products. (4) 
Reduced production of the product under consideration, in fact 
mirroring approach (2), not used in LCAs so far, but in inte¬ 
grated assessments on the effect of low-meat diets (Stehfest et 
al. 2009). (5) Land conversion emissions and delayed uptake 
due to land occupation (Muller-Wenk and Brandao 2010). 

For all approaches, the assumptions on reference case 
and time horizon are crucial (see Table 1). 

Approaches (1) and (3) assume that a certain starting value 
of production is not related to any land-use change, and can 
thus be maintained without conversion emissions. In ap¬ 
proach (3), historic conversion emissions due to an increase 
in production and land-use are spread out over the total 
production. Thereby, land-use change emissions are com¬ 
pletely dependent on the historic relative increase in produc¬ 
tion. This relation between production change and total 
production may change in the future, though. If the produc¬ 
tion in the start situation is large compared to the increase, 
conversion emission are on average very low. However, as 
assumed in approach (4), reducing the production would 
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Table 1 Comparison of different possible approaches and underlying assumption for including land-use and land-use change in LCAs 

Approach 

Description 

Assumed reference 

Time horizon 

Applied in 

(1) Continued 
current state 

Everything stays as it is now. 
This “approach” had not 
been applied explicitly, but 
implicitly this is why carbon 
implications of land-use 
have been ignored 
for so long 

Current state is assumed as 
reference. And as everything 
is assumed to stay in current 
state, no emissions or sinks 
need to be taken into 
account (obviously wrong) 

Irrelevant 

Most conventional LCAs until 
recently 

(2) Additional 

Assess the impact if additional 

Current state, or a baseline 

Very relevant, typically 

Most biofuel studies (e.g. 

production 

products would be produced 
compared to current state, 
or above a baseline 
projection 

development is assumed as 
reference. The choice of the 
baseline has limited 
influence on the result. 

conversion emissions are 
annualized over 20 or 

30 years. If longer time 
horizon is applied, 
conversion emissions 
decrease. 

(Searchinger et al. 2008)) 

(3) Historic land- 

Historic conversion emissions 

Situation at the beginning of 

Relevant, conversion 

Several LCAs for livestock 

use changes 

are attributed to the sector, 
and assigned to its total 
production. 

the historic period is used as 
reference. If production and 
land-use would not have 
increased, emissions would 
be zero. 

emissions depend on 
increase in production 
during the historic 
period. 

(Gerber et al. 2010; 

Ponsioen and Blonk 2011) 

(4) Reduced 

Assess the impact if less 

Current state, or a baseline 

Very relevant. Negative 

Method presented here 

production 

products would be produced 
compared to current state, or 
compared to a baseline 
projection, (can also be 
understood as “Missed 
Potential Carbon Sink” 
of land occupation) 

development is assumed as 
reference. The choice of the 
baseline has limited 
influence on the result. 

emissions due to avoided 
emissions and uptake due to 
land abandonment. Can be 

20 or 30 years (following 
(2)), or 100 years (see text 
for details). 

Implicitly applied in (Nguyen 
et al. 2010). 

(5) Transformation 

Separate treatment of 

Direct restoration of the 

Restoration time of natural 

Method is presented in 

and delayed 

transformation and land 

natural vegetation is 

vegetation (max 238 years) 

(Muller-Wenk and Brandao 

restoration 

occupation. Land 
occupation during one year 
results in a one-year delay 
of carbon uptake 
by regrowing natural 
vegetation 

assumed as reference, if 
land occupation would 
not be continued 

and 500 year cutoff time for 
C0 2 emission pulse. 

2010), no example 
application yet. No explicit 
statement on how to 
combine transformation and 
occupation. 


result in a decrease in land-use, a subsequent regrowth of 
natural vegetation at this location (e.g., forests, savannah 
tundra), and an associated uptake of C0 2 in this vegetation 
(Stehfest et al. 2009). In other words, a reduction in agricul¬ 
tural area leads to a larger carbon sink, and therefore one can 
argue that agricultural products are thus related to a “missed 
negative emission”. Using land for a certain product brings 
about that it cannot fulfill its potential as a carbon sink to 
mitigate GHG concentrations in the atmosphere. In this paper, 
we further elaborate approach (4), which has so far not been 
applied systematically in life cycle assessments. However, a 
recently published assessment on beef production systems in 
the EU in principle followed the concept, though assuming 
that the missed carbon sink would be identical to the emis¬ 
sions occurring during conversion from natural area to agri¬ 
cultural land, annualized over 20 years (Nguyen et al. 2010). 

All approaches (2)—(4) combine transformation and oc¬ 
cupation of land to one annualized indicator, and consider 
C0 2 fluxes as the primary indicator. Muller-Wenk and 


Brandao (2010) (approach 5) however, keep transformation 
and occupation separated. The C0 2 implications of land 
transformation for agricultural use in a certain year is de¬ 
scribed as C0 2 emissions from conversion, followed by a 
C0 2 uptake of the natural vegetation, which starts to regrow 
directly in the following year. Thereby an average stay of 
additional C0 2 from land transformation is computed. If 
land-use continues, the land occupation is described as a 
delay in the carbon uptake, resulting in one year prolonged 
stay of C0 2 in the atmosphere. Additionally, the authors 
compare the average stay of land-use related C0 2 to the 
average stay of fossil C0 2 . 

The aim of this paper is to elaborate approach (4), the 
methodology to include C0 2 implications of land occupa¬ 
tion (the “missed potential carbon sink”) in LCAs, and to 
apply it in a life cycle assessment of livestock products to 
derive total GHG emissions per kg product. We discuss our 
results and compare our approach to recent literature on the 
subject. Only by including the C0 2 emissions of land-use in 
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LCAs, the GHG effects of different products can be com¬ 
pared and be used to inform consumers on the consequences 
of their choices. 

2 Methodology 

The total GHG effect of a product is calculated as the sum of 
the emissions along the product chain according to conven¬ 
tional LCA (not including direct emissions from land-use 
change) plus the C0 2 emission or missed potential carbon 
uptake due to land-use occupation. 

GHGtotal = GHGLCAconventional T GHGmissedPotentialCarbonSi nk 

For the direct emissions along the production chain (GHG LCA 
conventional)? we use data from (Blonk et al. 2008), as they 
provide a consistent inventory for several livestock products, 
and also report area requirements per region (e.g., also for beef 
produced in the Netherlands they provide associated area 
requirements, e.g. for feed production in other regions, see 
Table 2). These data do not include any direct or indirect 
emissions due to land-use change or land occupation. 

The method for calculating the missed potential carbon 
sink is as follows: 

GHG ' m i sse dPot en ti a lCarbonSink 

L,R 

= 1/timehorizon x Area/ ?r x (CarbonSink/^) 

/=l,r=l 

whereby Area/ r is the agricultural area [in square meters] of 
land-use / (crop or grassland) in region r, needed per unit of 
product [in square meters per kilogram product], and Car- 
bonSink i r t is the carbon sink [in kilograms C0 2 per square 
meter] that occurs in region r when land-use / (crop or grass) 
is regrowing to natural vegetation (e.g., forests, or tundra) 
during t years. The time horizon is the time over which the 
potential C0 2 uptake is annualized. 


The “region” can be any geographic unit involved in the 
production process with characteristic current carbon con¬ 
tent and potential carbon sink, from small grid cells to world 
regions. This differentiation is very relevant, as the carbon 
stocks and thus potential sinks differ significantly across 
ecosystems. Likewise, the potential carbon sink also 
depends on the current land-use system of the product we 
are assessing, with grasslands often containing already more 
carbon than cropland systems. Additionally, the time t dur¬ 
ing which the C0 2 uptake is accumulated has to be defined. 
The carbon fixation is higher in the initial phases when the 
trees start to grow again, and gets smaller when the forests 
approach maturity. Finally, also a time horizon for allocating 
the carbon uptake needs to be defined. For the time horizon 
and the time t during which carbon is accumulated the same 
value should be applied, in the remainder we only talk about 
time horizon. For the process of carbon uptake, a period of 
100 years would be adequate, as by then the vegetation is 
coming close to its equilibrium state (e.g., Mila i Canals et 
al. (2007) mention a relaxation time of approx. 100 years). 
For biofuel studies, where mostly approach (2) is applied, 
which mirrors approach (4), the time horizon used for allo¬ 
cating emissions from land-use change is normally set to 
30 years (Searchinger et al. 2008), which is slightly longer, 
and thus a more conservative estimate than the 20 years 
suggested by IPCC for soil carbon processes (IPCC 1996). 
The 30 years period is also chosen as global GHG emissions 
have to be reduced strongly to achieve climate stabilization, 
and as reduction in the coming 30 years is both very difficult 
to achieve and crucial to avoid irreversible adverse effects 
from climate change (IPCC 2007). Here, we therefore 
explore results for both 30-year time horizon and 100- 
year time horizon, and also provide intermediate infor¬ 
mation for a 50 year time horizon. 

The area requirements per product in this paper are de¬ 
rived from LCA studies (Blonk et al. 2008), whereby dif¬ 
ferent production methods and locations are taken into 


Table 2 Area requirement of different livestock products in m 2 /kg product (Blonk et al. 2008) 


Livestock product 

Pasture 

Cropland 

Total 

Total area in SE-Asia 
and S-Am. 

Total area 
in N-Am. 

Total area 
in S-Am. 

Total area 
in Eur. 

Total area 
in Asia 

Total area 
in Austr. 

Beef, Brazil 

420.2 

0.0 

420.2 

420.2 

0.0 

420.2 

0.0 

0.0 

0.0 

Beef, Ireland 

54.6 

5.7 

60.3 

1.3 

1.5 

0.9 

56.0 

0.4 

1.5 

Beef, Netherlands (beef cattle) 

1.4 

13.3 

14.7 

1.4 

1.7 

1.0 

10.2 

0.4 

1.3 

Beef, Netherlands (dairy cattle) 

4.7 

2.6 

7.3 

0.6 

0.0 

0.6 

6.7 

0.0 

0.0 

Lamb, Netherlands 

30.3 

2.8 

33.1 

0.9 

0.0 

0.9 

32.2 

0.0 

0.0 

Calf, white, Netherlands 

0.9 

3.3 

4.2 

0.4 

0.3 

0.3 

3.5 

0.1 

0.0 

Pork, Netherlands 

0.0 

7.7 

7.7 

4.4 

0.3 

3.0 

2.6 

1.4 

0.0 

Chicken, Brazil 

0.0 

7.3 

7.3 

7.0 

0.3 

7.0 

0.0 

0.0 

0.0 

Chicken, Netherlands 

0.0 

4.6 

4.6 

2.6 

0.4 

2.6 

1.5 

0.0 

0.0 

Milk, Netherlands 

0.59 

0.32 

0.91 

0.07 

0.07 

0.05 

0.76 

0.02 

0.0 
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account. While agricultural products require the use of ag¬ 
ricultural area, the area requirements for other economic 
sectors will be much smaller. In this paper, we focus on 
agricultural products only, but the general principle, the 
calculation methods presented in this paper, can be applied 
to any product in a LCA calculation. 

The potential carbon sink is calculated using the IMAGE 
model (MNP 2006), per world region. IMAGE is an inte¬ 
grated assessment model to study global environmental 
change. Agricultural demand is calculated for 24 world 
regions, 5 livestock and 7 crop categories, and has been 
calibrated historically to FAO data. All physical land pro¬ 
cesses (like carbon cycling, crop production, animal graz¬ 
ing) are represented on a 30-min grid scale, and represented 
by process-based models with a monthly or annual time step 
(MNP 2006). land-use is allocated on the 30-min grid within 
a region until the actual or predicted production in that 
region is fulfilled. If less land is needed than in the time 
step before, some agricultural land is abandoned, and the 
potential natural vegetation can regrow there (for more 
detail see (Stehfest et al. 2009)). Carbon cycle processes in 
IMAGE are calculated by a modified version of the BIOME 
model (Van Minnen et al. 2000). Effects of land-use change 
on the carbon cycle are modeled dynamically; natural dis¬ 
turbances like fires are not modeled explicitly, but assumed 
to be captured in the equilibrium state of ecosystems. We 
carried out 30-, 50-, and 100 year-simulations where all 
parameters are kept constant, but reducing the production 
of beef, dairy, pork, sheep and goat, or poultry, respectively. 
In order to calculate the average effects for these five sec¬ 
tors, the experiment reduces the respective production to 
zero. Thereby, we derive the decrease in agricultural area, 
the increase in carbon stock, and thus the average carbon 
sink per area for the different livestock products. By follow¬ 
ing this approach, we calculate the consequences of an 
average product. However, the first unit of product not 
consumed any more would have different consequences 
than a later unit not consumed any more, depending, e.g., 
on whether intensive or extensive cattle systems would be 
abandoned first. 

3 Results 

The area requirements of different livestock products in dif¬ 
ferent world regions as derived from (Blonk et al. 2008) are 
shown in Table 2. As described in the Section 2, the IMAGE 
model has been used to derive the potential carbon sink for 
different time horizons and animal products. This potential 
carbon sink per world region, corresponding to the regions in 
Table 2 and in Blonk et al. (2008), is shown in Table 3. 

Combining these areas and the potential carbon sink per 
square meter from the IMAGE experiments results in the 


potential carbon sink per livestock product (see Section 2). 
In Table 4, we show the direct GHG emissions from the 
production chain (derived from Blonk et al. (2008)), the 
missed potential C0 2 sink for a 30- and 100-year time 
horizon, and their sum. 

For grazing animals, applying a time horizon of 30 years, 
the missed potential C0 2 sink ranges from 2 to more than 
200 kg C0 2 /kg product, 1 and exceeds the conventional LCA 
GHG emissions by far. For pork and poultry, applying a 
time horizon of 30 years, the numbers are in the same order 
of magnitude, which means that in these cases, the LCA 
totals including C0 2 emissions from land-use lead to LCA 
results that are about twice the results of conventional LCA. 

Using a 100-year time horizon for the missed potential 
carbon sinks results generally in smaller results than when 
applying a 30-year time horizon. This is not surprising as the 
carbon sink of regrowing forests on abandoned areas is 
declining after the trees have finished their initial strong 
growth in the first decades. Nevertheless, the missed poten¬ 
tial carbon sink results for a 100-year time horizon are still 
highly significant. 

As a comparison, we also calculated the missed potential 
carbon sink for protein rich, plant-based products (Table 5). 
The LCA totals including C0 2 related to land-use vary be¬ 
tween 2.4 kg C0 2 -eq/kg for Tempeh, and 3.78 kg C0 2 -eq/kg 
for Tofu when applying a 30-year time horizon. For the 
potential carbon sink for a region per square meter, we use 
the average values taken from Table 3 here. All results for 
livestock and plant products are shown in Fig. 1. 

In the above calculations the IMAGE model was only 
used to determine the potential carbon sink per area, for 
different regions, which were then combined with the area 
requirements from an LCA study (Blonk et al. 2008). We 
also used the entire IMAGE model output to calculate the 
potential carbon sink per unit of product (Table 6). These are 
world average numbers, and therefore they differ from the 
results in Table 4 for some livestock products. For beef, they 
are in the wide range of the results for different systems, 
for other products the world average numbers are signif¬ 
icantly higher, as they include much less intensive pro¬ 
duction systems across the world. 

In order to make our approach applicable in LCA studies, 
we included information on the development of the missed 
potential sink and the increase in the carbon stock (see 
Table 3). 


1 One kilogram of C0 2 corresponds roughly to the emissions from a 
car driven for 5 km, assuming an average European car emitting 160 g 
C0 2 ,/km from fuel combustion (VCD—Verkehrsclub Deutschland 
2008; BFE Bundesamt fur Energie 2009) plus 40 g C0 2 emissions 
for vehicle and fuel production (Lane 2006). Applying these calcula¬ 
tions to the LCA results here, 1 kg of Brazilian beef would translate to 
driving an average European car over 1,600 km, 1 kg of Dutch chicken 
to 31 km (using the results using the 30-year time horizon), and so on. 
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Table 3 Average potential car¬ 
bon sink per region and product 
type as calculated by the 
IMAGE model (above). 
Averages across these scenarios 
for the potential carbon sink, and 
the additional carbon in the 
vegetation below 


Regional breakup was selected 
to fit regions in Blonk et al. 
(2008) 


Table 4 GHG emissions from 
conventional LCA, missed 
potential carbon sink and total 
GHG emissions for different 
livestock products and origins 
[kg C0 2 -eq/kg product] 


GHG LCA data from (Blonk et 
al. 2008). Time horizon is 30 and 
100 years, respectively 




World 

Avg 

N- 

America 

Latin 

America 

Europe 

Asia 

Australia 

Time horizon 



kg C0 2 /nr 

’/year 



Less beef 

30 

0.31 

0.42 

0.66 

0.47 

0.16 

0.01 


50 

0.29 

0.42 

0.58 

0.47 

0.17 

0.00 


100 

0.20 

0.33 

0.40 

0.38 

0.13 

0.00 

Less dairy 

30 

0.51 

0.49 

0.85 

0.51 

0.35 

0.64 


50 

0.46 

0.50 

0.68 

0.53 

0.36 

0.58 


100 

0.33 

0.38 

0.45 

0.43 

0.28 

0.43 

Less pork 

30 

0.66 

0.45 

1.07 

0.53 

0.58 

0.95 


50 

0.56 

0.47 

0.77 

0.54 

0.55 

0.76 


100 

0.40 

0.38 

0.51 

0.44 

0.40 

0.55 

Less sheep and goat 

30 

0.53 

0.94 

1.14 

0.57 

0.36 

0.64 


50 

0.46 

1.16 

0.75 

0.56 

0.37 

0.58 


100 

0.32 

0.25 

0.39 

0.45 

0.30 

0.43 

Less poultry 

30 

0.69 

0.42 

1.01 

0.56 

0.67 

0.96 


50 

0.58 

0.46 

0.71 

0.55 

0.60 

0.71 


100 

0.41 

0.35 

0.47 

0.47 

0.44 

0.49 

Average across these scenarios 












kg C0 2 /m" 

’/year 



Potential carbon sink 

30 

0.54 

0.54 

0.95 

0.53 

0.42 

0.64 


50 

0.47 

0.60 

0.70 

0.53 

0.41 

0.52 


100 

0.33 

0.34 

0.44 

0.43 

0.31 

0.38 





kg C/m 

2 



Additional carbon in 

30 

4.44 

4.44 

7.74 

4.31 

3.48 

5.23 

vegetation 

50 

6.41 

8.20 

9.51 

7.23 

5.56 

7.15 


100 

9.09 

9.23 

12.12 

11.82 

8.47 

10.40 


Product Product GHG Missed Total GHG Missed Total GHG 

Region LCAconv potential emissions potential emissions 

carbon sink carbon sink 

30-year time horizon 100-year time horizon 


kg C0 2 -eq/kg 


Beef 

Brazil 

59.0 

276.1 

335.1 

166.4 

225.4 

Beef 

Ireland 

38.3 

27.4 

65.7 

22.1 

60.4 

Beef 

4.9 

Beef cattle 
20.8 


Netherlands 

15.9 

6.2 

22.1 

Beef 

2.8 

Dairy cattle 
11.7 


Netherlands 

8.9 

3.5 

12.4 

Lamb 

Netherlands 

16.3 

19.3 

35.6 

14.8 

31.1 

Calf 

Netherlands 

6.3 

2.0 

8.3 

1.6 

7.9 

Pork 

Netherlands 

4.5 

5.6 

10.1 

3.4 

7.9 

Chicken 

Brazil 

2.6 

7.2 

9.8 

3.4 

6.0 

Chicken 

Netherlands 

2.6 

3.6 

6.2 

2.1 

4.7 

Milk 

Netherlands 

1.2 

0.5 

1.7 

0.4 

1.6 
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Table 5 Area requirement, conventional GHG LCA data, missed potential carbon sink, and total LCA GHG effect for the production of different 
meat-alternatives, using a 30-year time horizon 


Product and land 
of production 

Cropland 

Total area 
N-Am. 

Total area 
S-Am. 

Total area 
Eur. 

Total area 
Asia 

GHG 

LCAconv 

Missed potential 
carbon sink 

Total 

GHG emissions 





m 2 /kg 



kg C0 2 -eq/kg 


Tofu, Netherlands and Belgium 

3.036 

1.821 

0.455 

0.455 

0.304 

2.0 

1.8 

3.8 

Tempeh, Netherlands 

2.260 

1.356 

0.339 

0.339 

0.226 

1.1 

1.3 

2.4 

Quorn with 4 % egg white, UK 

1.101 

0.044 

0.420 

0.634 

0.003 

2.6 

0.8 

3.4 


Areas and GHG LCA data from Blonk et al. (2008) 


4 Discussion and conclusions 

Up to now, life cycle inventories for GHG emissions of products 
have largely ignored area demands and the emissions and po¬ 
tential sinks related to land occupation. As a consequence, LCA 


results do mostly depend on process chain GHG emissions. 
However, whether agricultural production requires a lot or little 
land has major implications for GHG emissions and climate 
change, revealing a major loophole in hitherto existing LCAs, 
especially for agricultural products. With the recently increased 


Fig. 1 GHG emissions for the 
various products calculated in 
Tables 4 and 5, using a 30-year 
time horizon. The chart shows 
the GHG emissions according 
to conventional LCA, the 
missed potential carbon sink, 
and the resulting total GHG 
effect 
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Table 6 Global average missed potential carbon sink for different livestock products and time horizons [kg C0 2 -eq/kg product] calculated entirely 
based on IMAGE results, i.e., IMAGE area requirements, not using numbers from Blonk et al. (2008) 


Time horizon 

Beef 

Milk 

Pork 

Poultry 

Sheep and goat 

30 

192.1 

4.9 

18.0 

15.9 

141.9 

50 

174.7 

4.3 

14.9 

13.0 

122.3 

100 

124.0 

3.1 

10.4 

9.1 

85.5 


attention for the climate effects of especially livestock produc¬ 
tion, and increased demand for environmental labeling, it is 
necessary to develop methodologies for GHG emissions related 
to land occupation, and to include them in LCAs. 

Although its consequences for emissions and potential 
sinks have largely been ignored, some LCA studies and similar 
approaches address the occupation of land as an individual 
indicator. The ecological footprint method does account for 
that in global hectares, or virtual “earths” (see, e.g., Kitzes 
et al. (2008) and Ewing et al. (2008)). Some life cycle 
impact assessments address several impacts of land occu¬ 
pation, from nutrient and water cycling to biodiversity 
(Mila i Canals et al. 2007). 

And recently, Muller-Wenk and Brandao (2010) have 
presented an approach to include carbon transfers between 
land and the atmosphere into LCA. We therefore discuss our 
study in relation to their methodology. First of all, they 
explicitly distinguish between the transformation and the 
occupation of land, while our methodology, and also 
approaches (2) and (3) combine the transformation and 
occupation into one annualized indicator. While the distinc¬ 
tion is of course accurate from a conceptual point of view, it 
does not reflect the intended applications of many LCAs, 
and may thus be misleading. If land transformation is only 
attributed to the first year’s harvest of, e.g., soy beans, while 
all later harvests are only linked to land occupation, i.e. 
postponed relaxation, this will lead to very different impact 
indicators for otherwise identical products (the soy beans). 
Therefore, we claim that the transformation emissions need 
to be distributed over the occupation period in some way or 
another. Muller-Wenk and Brandao (2010) do not offer a 
methodology for this, but in principle one could use their 
numbers and add up a certain fraction of transformation 
(e.g., 1/30) to every year of occupation. The second major 
difference is that Muller-Wenk and Brandao (2010) do ac¬ 
count for the average stay of C0 2 in the atmosphere: They 
calculated 157 years for fossil C0 2 , and 0.5 times the time 
needed for the natural vegetation to regrow (“relaxation 
time”) for emissions due to land-use change. They then 
calculate the fossil fuel equivalent C0 2 emitted to the air, 
by relating these two different times to a duration factor. The 
underlying reasoning is that any carbon emitted by land will 
be taken up again by the land once the occupation stops. In 
their argumentation, the end point of the alternative 


situations should be identical; instead of a “missed carbon 
sink”, they would talk about a “postponed carbon sink”. 
However, we argue that such an approach could again be 
misleading for the application of LCAs in the context of 
climate change mitigation. A molecule of C0 2 now has a 
distinct different effect on the climate system than a carbon 
molecule 30 years later (see O’Hare et al. (2009) for dis¬ 
cussion on time accounting), and for climate stabilization, 
the largest challenge is to reduce emissions and C0 2 -levels 
in the atmosphere in the coming few decades (IPCC 2007). 

In terms of carbon stocks associated to certain land trans¬ 
formations, our numbers are consistent with the sources 
cited by Muller-Wenk and Brandao (2010), as the IMAGE 
carbon cycle model has also been evaluated against global 
carbon stock inventories (Klein Goldewijk et al. 1994). For 
example, for Latin America, the carbon sink over a 100 year 
period ranges from 0.4 to 0.5 kg C0 2 /m /year (see Table 3), 
corresponding to a total stock of 40-50 kg C0 2 /m , or 110- 
136 t C/ha. This is just slightly lower than their value for 
“tropical forests” (151 t C/ha), which probably reflects the 
fact that not all of the land-use change in Latin America is 
occurring in tropical forests. 

As also discussed in Section 2, we calculated our results 
for a 30- and 100-year time horizon, consistent with current 
literature. While the choice of the time horizon remains 
subjective, reflecting preferences in time, most biofuel stud¬ 
ies apply a 30-year time horizon (Searchinger et al. 2008), 
and many other studies use multiple time horizons, often 
including 100 years, as also global warming potentials are 
often applied over 100 years. Therefore, we decided to 
present our results for a 30- and 100-year time horizon. As 
we distribute the carbon sink equally over all years, and as 
the carbon uptake of regrowing vegetation is leveling off, 
the annual missed potential carbon sink is decreasing 
with an increasing time horizon. Although we are aware 
of recent efforts to include temporal aspects in LCAs 
(Kendall et al. 2009; Levasseur et al. 2010; O’Hare et 
al. 2009), we do not include such methods here, as they 
are less relevant for the (continuous) carbon uptake than 
for initial emission peaks. We also do not include a 
discounting of emissions, following the USEPA regula¬ 
tions on biofuels (USEPA 2010). 

As we estimate the GHG implications of land occupation 
by calculating the missed potential carbon sink, analyzing 
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the uptake of C0 2 if the land was no longer used for 
agricultural production, our approach might be considered 
as a consequential LCA approach. While attributional LCAs 
use average data on a given situation of production, conse¬ 
quential LCAs analyze the effects of (small) changes in the 
production of goods. 

Additionally, it might be interesting to discuss the rela¬ 
tion of our approach to assessments of NPP on natural and 
agricultural land. Weidema and Lindeijer (2001) use the 
NPP (net primary production) differences of natural vegeta¬ 
tion and occupied land for the calculations of occupation 
impacts during the production of agricultural goods, and 
also the concept of HANPP (human appropriation of net 
primary production) assesses the difference in NPP due to 
human influence (Haberl et al. 2007). While the assessment of 
NPP is suited to assess the general impact of human land-use 
on ecosystem productivity, “biotic resources” and “life-sup- 
port functions of natural systems” (Weidema and Lindeijer 
2001), it is, however, not suited to assess the GHG and climate 
implications of land occupation (most of the NPP does not 
constitute a net exchange of the entire ecosystem with the 
atmosphere, but the C0 2 fixed via NPP is mostly respired 
again via heterotrophic respiration). 

The novel approach presented here to include the C0 2 
implications of land occupation via the missed potential 
carbon sink provides the possibility to account for this 
important contribution to the overall GHG balance. Al¬ 
though the process of carbon uptake is temporary, and stops 
when the new equilibrium of the regrowing natural vegeta¬ 
tion is reached, this regrowth might be especially relevant 
under ambitious climate policy, where strong emission 
reductions are needed during the next few decades. 

Being simple by design, there are several uncertainties 
related to this approach: We calculated the average missed 
potential carbon sink for all current consumption. However, 
land-use systems are very heterogeneous, and thus the first 
unit of product not consumed any more is likely to have 
different consequences than a later unit not consumed any 
more. For example, the least intensive and most land con¬ 
suming production systems may be abandoned first. The 
same heterogeneity applies for the natural vegetation that 
would come in place of the abandoned land. Depending on 
the location, the carbon content of the regrowing vegetation, 
and thus the potential sink varies a lot, with typical least 
intensive systems on marginal grasslands leading to a 
smaller carbon sink. However, as the dynamics of decreas¬ 
ing production and its location are very uncertain, it was 
decided to follow this average approach. 

Furthermore, we assumed that the entire land currently 
occupied for a certain product will be converted into natural 
vegetation. Doing so keeps the method simple and easily 
applicable, and follows the assumption of a clear separation 
between the environmental and the economic system in 


LCAs. However, the consequences of reduced production 
may feed back into the economic system, and via price 
effects the abandoned area may in fact be smaller. Some 
developments for consequential LCAs try to capture these 
feedback by applying global agricultural models, e.g. for 
expanding production of biofuels (Kloverpris et al. 2010; 
Kloverpris et al. 2008; Taheripour et al. 2010). In future, more 
detailed analysis, these feedbacks might be taken into account, 
but they are also highly uncertain (Stehfest et al. forthcoming). 

If livestock grazing or cropland for feed is removed, the 
regrowth of the vegetation naturally occurring at this location 
(e.g., forest, savannah, or tundra) is not inevitably the only 
possibility, and the scenario “natural vegetation replaces agri¬ 
cultural areas if possible” is somewhat arbitrary. There may be 
other alternative uses, like biofuels or photovoltaics, also 
leading to (possibly even stronger) benefits in the GHG bal¬ 
ance than the carbon uptake due to regrowth of vegetation. 

Finally, it has to be noted that some of the extensive 
grazing land is not only having a low potential carbon sink 
as described above, but that for some of these areas grazing 
is the only way it can contribute to food production, as it is 
not suitable for growing crops. Consequently, in certain 
areas replacing extensive, grazing cattle systems by more 
efficient, intensive, feed-based cattle systems is not a rea¬ 
sonable option, as it would increase the pressure on areas 
suitable for crop production elsewhere. Contrary, in regions 
where grasslands are also suitable for feed crop production, 
switching to more intensive systems will lead to less land 
occupation, though the replacement of grass by feed crops is 
limited by beef physiology. Compared to grazing ruminants, 
meat from monogastrics in general involves less land occu¬ 
pation, and most plant-based products have even markedly 
lower land-use implications. However, the impacts of soy- 
based products like tofu can be more than half as high as 
those of intensive chicken breeding, as yields of soybean are 
much lower than the yield of typical feed crops like maize. 

The new integrated LCA approach presented in this paper 
should lead to more realistic GHG emission results, as 
products consuming much land for their production will 
get additional C0 2 -emissions added to their LCA balance, 
representing the missed potential carbon sink that could be 
realized if the product would not be consumed. 
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